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Abstract. We investigate the clustering properties of dynamical Dark Energy even in 
association of a possible coupling between Dark Energy and Dark Matter. We find that 
within matter inhomogcneitics. Dark Energy migth form voids as well as overdensity 
depending on how its background energy density evolves. Consequently and contrarily 
to what expected. Dark Energy fluctuations are found to be slightly suppressed if 
a coupling with Dark Matter is permitted. When considering density contrasts and 
scales typical of superclusters, voids and supervoids, perturbations amplitudes range 
from l^^l ~ 0(10"^) to \5^\ ~ 0(10"^) indicating an almost homogeneous Dark Energy 
component. 
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1. Introduction 

Several observations made over tlie recent years, related to a large extension to Large 
Scale Structures (LSS) and anisotropies of the Cosmic Microwave Background (CMB) 
as well as the magnitude-redshift relation for type la Supernovae [1], have given us a 
convincing picture of the energy and matter density in the Universe. 

Baryonic matter accounts for no more than 30% of the mass in galaxy clusters while 
the existence of a large clustered component of Dark Matter (DM) seems now firmly 
established, although its nature is still unknown. However, they contribute to the total 
energy density of the Universe with only a few percent and about 25% respectively. 

No more than another few percent could be accounted for by massive neutrinos, 
but only in the most favorable, but unlikely case. According to [2j the total mass of 
neutrinos cannot exceed the limit of 1.43 eV (see, however, [3] for a recent analysis on 
neutrino mass limits in coupled dark energy models). A very small part (10~^) of the 
total energy density is due to massless neutrinos and CMB radiation. 

The model suggested by observations is only viable if the remaining 75% is ascribed 
to the so-called Dark Energy (DE) responsible for the present day cosmic acceleration. 

Although strongly indicated by the observations, the existence of DE is even 
more puzzling than DM. It can be identified with a cosmological constant or with a 
yet unknown dynamical component with negative pressure. On the other hand, its 
manifestation can be interpreted as a geometrical property of the gravity on large scales 
resulting from a failure of General Relativity (GR) on those scales (see |1] for a review). 

Within the context of GR, as an alternative to the cosmological constant, DE is 
usually described as a self-interacting scalar field or a cosmic fluid with negative pressure 
(see and references therein). It is usually assumed that density perturbations of DE 
play a negligible role in the structure formation because of its very small mass ~ H {H 
being the Hubble parameter). Accordingly, perturbations should appear only on very 
large scales (>100 Mpc) and are bound to be linear so that rates of structure formation 
and their growth are influenced by DE only through the overall cosmic expansion [6]. 

Nevertheless, this assertion remains questionable and its validity has been subject 
to recent debate. Then, the key question is whether DE actively participates in 
the clustering and virialization processes possibly developing non-linearity on relevant 
scales. Some attempts to solve the problem have been done in [71 [HI O [10], sometimes 
adopting a phenomenological approach parametrizing the observables associated with 
DE clustering or in the context of coupled DE and other non-trivial DE models where 
clustering is expected to be more probable. 

An intriguing result was obtained by Dutta & Maor [9]. They numerically 
studied the clustering of DE within matter overdensities showing that DE tends 
to form underdensities or voids in response to the gravitationally collapsing matter 
(similarly, DE overdensities are expected in correspondence of matter underdensities). 
On supercluster scales they found ~ (9(10^^) when 6m ~ C'(l) which could 
be relevant to observations {6m and 6^ being the density contrast of matter an DE 
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respectively). However, their analysis was limited to linear spherical perturbations and 
the simplest class of DE models in which DE is ascribed to a ligth scalar field (j) slowly- 
rolling down its pontential V^(0), minimally coupled to matter and gravity. 

A similar problem was investigated by Mota et al [10] with a different approach 
providing an analytical approximation for 6^ valid both in linear and nonlinear regime 
{6m > in uncoupled scalar field DE models. However, in spite of the qualitative 

agreement of their results with those of Dutta & Maor [9], they found 15^1 ~ O(10~^) 
when the same linear scales are considered. 

In the present paper we extend the analysis of Dutta & Maor [9] taking into 
account a possible coupling between DM and DE and considering a wider class of DE 
tracking potentials. Then, we numerically study DE clustering on those scales for which 
Sm ~ C'(l) today, i.e. supercluster, void and supervoid scales, and show that, coupled 
as well as uncoupled DE develops inhomogeneities which amplitudes are consistent with 
the findings of Mota et al [10]. Formation of DE voids (overdensities) is however related 
to how the background energy density of DE evolves and not only to the presence of 
matter overdensities (underdensities) as claimed in the previous works. It is also shown 
that, as expected, the growth of matter fiuctuations is suppressed as soon as DE starts 
to drive the cosmic acceleration while 5^ — as DE, in the models here considered, 
asymptotically approaches a cosmological constant. 

The plan of the paper is as follows: in Section [2] we describe our model and give 
the linearized equations for matter, DE and metric perturbations which are derived in 



Appendix A Numerical results concerning perturbation evolution are presented and 



discussed in Section |3l We summarize and conclude in Section HI 



2. The model 



The essential feature of a scalar field (f), in order to yield DE and drive the cosmic 
acceleration, is its self-interaction through a potential V{(j)). 

In addition to self-interaction, a scalar field can in principle be coupled to any 
other field present in nature. Couplings to ordinary particles are strongly constrained 
by observational limits on violations of the equivalence principle but limits on the DM 
coupling are looser (constraints on coupling for specific models were obtained in [TTl [T2] 
from CMB, N-body simulations and matter power spectrum analysis). If present, DM 
coupling could have a relevant role in the cosmological evolution affecting not only the 
overall cosmic expansion but also modifying the DM particles dynamics with relevant 
consequences for the growth of the density perturbations in both the linear and nonlinear 
regime (e.g., on halo density profiles, mass function and its evolution) [TT| [T3llT^ . Here 
we consider one of the most popular models where a coupling between DM and DE is 
present, namely coupled DE [13] (for different DE-matter interactions see [T5]). 

According to general covariance, the sum of the individual stress energy tensors 
~ ^' ^) must be conserved so we can write: 



Voids and overdensities of coupled Dark Energy 



4 



Here is the co variant derivative, T(^dm) indicates the trace of the DM stress energy 
tensor and C is a constant which parametrizes the strength of the DM-DE interaction. 

In a spacetime described by a metric g^u with signature ( — h ++), the DE stress 
energy tensor takes the form: 

Ti;^, = d^ct^dA - 6i: i^ld^<jyd''<f> + v'^ (2) 

while baryons and DM are well described as non-relativistic pressureless perfect fluids 
with: 

= PiU^i^u^ii) (3) 
where pi is the energy density of the component i = b, dm and u'^--^ is its 4-velocity. 

As we are interested in spherical perturbations around a spatially fiat Friedmann- 
Robertson-Walker (FWR) background, only spherically symmetric spacetimes are 
considered for which the most general line element in comoving coordinates is: 
ds^ = -dt^ + U{t, r)dr^ + V{t, r) (d^^ ^ gj^2 ^^^2^ ^ 

where U{t,r) and V(t,r) are general functions [T6] . 



2.1. Background and perturbation equations 

As usual, equations for background and perturbation evolution follow from and 
Einstein's equations. In this Section we only state the full set of equations while their 
derivation is detailed in Appendix A Working in the synchronous gauge, we redefine 
the metric functions as follows: 

W(t,r) = a(t)V«*''-) 

V(t, r) = rVt) V^^*'") (5) 

Here a{t) is the scale factor of the homogeneous background while ( and ip represent 
deviations from homogeneity. In the following, however, metric perturbations will be 
described by x = C + 2^'? the only combination which is relevant to the equations of 



motion (see Appendix A for more details). We also decompose 0, pdm and pb as the 



sum of an unperturbed part, denoted by , and a perturbed one: 

<P{t,r) = ^{t) + 6<j){t,r) 

Pdm{t, r) = Pdm(t) + Spdmit, t) 

Pb{t, r) = pbit) + 6pb{t, r) (6) 

and, for each component i, we define the density parameter f2j = pi/pcr, the density 
contrast 6i = Spi/pi and 9i = div Vj {pcr and Vj being the critical energy density 
and the coordinate velocity respectively). We assume perturbations to be small so 
that linear approximation applies. Further, only radial motions are considered, i.e. 
|v| = v'"" = dr/dr. 

Hereafter we make use of conformal time r related to cosmic time t via the equation 
dt = adr. Derivatives with respect to r are denoted with an overdot and the conformal 
Hubble function Ti = a/ a = aH. 
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From ([2]), it is possible to work out the background energy density of DE, its 
pressure P^, the corresponding perturbations as well as the radial velocity: 



2a2 
dJ 



+ V'l 



V'i 



(7) 



where V = V{(f)) and ' denotes the derivative with respect to (p. 

With the above notation and definitions, the equations for background evolution 
read: 

37^2 ^ ^p^^ + pb + p^] 
I + 27^0 + aW = Cpdma'^ 

Pdm + 'i'^Pdm = —Cpdm<t> 

k + snpb = (8) 
while linear perturbations evolve according to: 

6(f) + 2m(j) - V^6(f) + aW'6(f) + x0 = Cddmpdma^ 

5dm + Odm + X= -C6(j) 
0dm + ^Odm = C ((j)6dm + 

6b + eb + x = o 
9b + n9b = o 



X 



A 2 








(9) 



Note that, in the absence of coupling, C = 0, DM particles and baryons follow the same 
dynamics and can be used to define the synchronous coordinates and therefore have 
zero peculiar velocity {9dm = 9b = 0). 9i's equations are not needed anymore and the 
set of equations ([8]) and Qj then reduces to that given in Dutta & Maor [9]. 

On the other hand, as widely discussed in [13], coupling modifies the dynamics of 
DM particles and an equation for 9dm is therefore needed: although 9dm = initially, it 
can not remain null for all times because of the term C (\/'^6(f) + (f)9dm^ ■ As a consequence 
baryons and DM develop a bias b, i.e. Sb = bSdm- Anyway, it is still possible to eliminate 
the variable 9b reducing the number of equations by one. 

It is also worth mentioning that, unlike the uncoupled case, in the presence of 
coupling. Universe goes through an evolutionary phase named 0-matter dominated era 
((/)MDE) just after matter-radiation equivalence. In this period the scalar field behaves 
as stiff matter {P^/ p^ = 1) having a non-negligible kinetic energy which dominates over 
the potential one. After this stage, the usual matter era follows before entering in the 
accelerated regime. 
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In order to perform a numerical integration, it is most convenient to work in the 
Fourier space. Setting 9^ = 0, equations ([9]) then become: 

6(f)k + 2H6(j)k + k'^S(j)k + a^V"5(j)k + Xk4> = C5dm,kPdma^ 

Sdm,k + Odm,k + Xk = —C5(j)k 

(^dm,k + 'H6dm,k = C ((j)6dm,k — k'^^4>kj 

kk + Xk = (10) 



Xk + nxk - -H' 



^dmSdm,k + ^bSb,k + ^(f) 5^,k 



<f>,k 







where the index k denotes the Fourier-components with wavenumber k. 

After having numerically evolved the above equations, their solutions are then 
Fourier transformed back to the coordinate space. 



2.2. Potential 

The present analysis is based on the assumption that DE is a self-interacting scalar field 
0. Two self-interaction potentials are considered [T7] : 

^(0)=A°+70" RP (11) 

or 

1/(0) = (A°+70")exp(47r^) SUGRA (12) 

ml 

{nip = is the Planck mass) admitting tracker solutions and initially introduced to 

ease the fine tuning and coincidence problems. In the absence of DM-DE coupling, RP 
yields quite a slowly varying w{a) = Pci,/p^ state parameter. On the contrary, SUGRA 
yields a fastly varying w{a). Although coupling causes a w{a) behavior significantly 
different from the uncoupled case, one could again consider these potentials as examples 
of rapidly or slowly varying w{a). 

For any choice of the energy scale A and the positive parameter a, the above 
potentials yield a fixed Vt^. Here we prefer to use A and Vl^ as free parameters; the 
related a value is then suitably fixed. A values are chosen according to the constraints 
given in [18]. 

The use of the above potentials also permit to verify the generality of the Dutta 
& Maor [9] results obtained for different ones, namely the mass potential V = m?(fp' /2 
and the double exponential potential [19]. 

In both uncoupled and coupled cases, the above potentials admit a final de Sitter 
attractor {Vt^ = 1, = —1) toward which universe asymptotically evolves. The behavior 
of the scalar field in this last stage is however different in the two cases. In the SUGRA 
model, driven by damped oscillations, the field approaches the minimun of its potential 
rather than asymptotically approaching infinity as in the RP case. As we will see in the 
next section, perturbations reflect the same behaviors. 
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2.3. Initial conditions 

As in Dutta & Maor [9], perturbations in baryons and DM are initially taken to be 
gaussian: 

5drn{ri,r) = A{Ti)e ^ 6b{Ti,r) = b6dm{ri,r) (13) 

(here is some initial time in the matter era). The ratio between fluctuation amplitudes 
in baryon and DM is prescibed by linear theory |13j : 



37 ^dm + A(3^y/rik 



(14) 



where (3 = 3 / IQirnipC is the adimensional coupling parameter, /x = {6dm/ ^ dm) I {o./ a)., 
7 = 1 + 4/32/3 and Vtk = 0V2aVcr- 

The initial perturbation amplitude A{Ti) is chosen such that the mean value of the 
total matter density contrast 6 = {6pb + Spdm)/{pb + Pdm) within the comoving radius 
R = a at the present time tq is: 

J^ 5{r,To)r'^ dr , 
/(f r'^dr 

when dealing with supercluster and voids scales (~ 10 — 25h^^ Mpc). Slightly smaller 
\6m\ will be considered for supervoid scales (see next section). 

We assume the shape of matter perturbations to be only slightly affected during 
their evolution so that 6dm and 6b are still well approximated by a gaussian at tq. We 
will see in the next section that this assumption is conflrmed by numerical results. Thus, 
setting 6 = N6dm = NAe ^we have: 



6{r <R,t) = ^NA 



V^a fR\ _ 

erf — — e 

2 R KaJ 



where = [h Vlb + ildm)/i^b + ^dm)- For R = a and r = tq it follows A ~ 1.76/A^ 
when (5 ~ 1. 

Dutta & Maor [9] set the initial redshift to Zi = 35 assuming 6(j), 6(j) = 0. However, 
DE fluctuations would have evolved since the earlier stages of the Universe and, although 
small, they are expected to differ from zero at Zi being then in tracking regime. With 
their choice for the initial conditions, Dutta & Maor [9] flnd that, initially, DE has a 
weak tendency to collapse. However, this effect has to be understood as due to the 
fact that such initial values do not lie on the tracker solution of the perturbed scalar 
fleld equation. This is clearly illustrated in Figure [1] where the evolution of the density 
contrast at the center of a matter overdensity is plotted for different choices of the initial 
conditions. In the top panel we compare the behaviors of 5,^ when choosing 6(f), 6(f) = 
at different initial redshifts, i.e. Zi = 100, 500, 1000. Notice the initial tendency of DE 
to cluster. However, for different choices of 6(f> and 6(f> the initial behavior of 6^ will be 
completely different. This is shown in the bottom panel, where we permit initial values 
to differ from zero. An inspection of the plots, shows that independently from the 
initial conditions chosen, after an initial transient period, 5^ always settles on the same 
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Figure 1. Evolution of the DE density contrast, 6^, at the center of a spherical 
matter overdensity for different initial conditions. In the top panel we set (5</>, Scj) = 
at different initial redshifts (zi — 100, 500, 1000) while in the bottom panel we plot the 
behavior of when the initial values of (50 and Scj) differ from zero. Independently 
from the initial conditions chosen, after an initial transient period, always settles 
on the same (tracker) solution. 



(tracker) solution. Plots are given for a specific RP model but the same conclusions are 
reached when considering different cases. 

Based on the above considerations, we choose initial conditions such that 50 
arranges itself on the tracker solution already at z = 100 (e.g. 6(f), 6(j) = 0, Zi = 1000). 

Finally, we set xl'^j? = assuring that matter perturbations initially expand with 
the Hubble flow. 

3. Results 

In this section we present the results of numerical runs. In order to better understand 
the behavior of the perturbations around the present time, we stopped our runs at some 
time in the future. Results are then shown in the redshift range 100 < z < —0.99 and, 
for a = 20h~^ Mpc if not otherwise specifled. All models considered are spatially flat, 
have adimensional Hubble parameter h = 0.70, Qf, = 0.046 and Qdm = 0.234. 

Let us start flrst with the uncoupled case. Figure [2] compares the time evolution 
of the DE density contrast, 6^, at the center of matter overdensity (r = 0), in RP 
and SUGRA (uncoupled) cases. Both models present the same quahtative behavior 
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Figure 2. Evolution of the DE density contrast, 6^, at the center of a spherical 
matter overdensity in RP (dotted line) and SUGRA (solid line) uncoupled models. 
The comoving scale of the perturbation is tr = 20h~^ Mpc and the mean matter 
density contrast within the radius i? = cr is 5 = 1 at the present time. 



until DE starts to dominate the cosmic expansion and approaches the Plack mass rrip. 
Differences in the late time behaviors are then to ascribe to the exponential term in (fT2ll 
which dictates the late evolution of (p in SUGRA model as explained in the previous 
section. 

Just as Dutta & Maor [9] and Mota et al we find that DE tends to form voids in 
correspondence of matter overdentisities obtaining |5<^| ~ O{10~^ — 10~^) in late matter 
era. While this is consistent with the findings of Mota et al [1^, it does not agree 
with those reported by Dutta & Maor [9], i.e. \6^\ ~ O{10~^). Following the same 
arguments by Dutta & Maor [9], voids formation is understood as a 'drag effect' due to 
the slower expansion of the regions with matter overdensities. Matter collapse lowers 
the local value of H, reducing the Hubble damping to the scalar field. Therefore, in 
those regions, rolls down its potential sligthly faster increasing its background value 
and by the quantities 50 and 50 respectively. Noticing that 0, 0, 50, 50 > while 
V' < 0, a local void in DE, i.e. 5p^ = 6pv + Spk < 0, will form when the local variation 
in the potential energy density, 6pv = V'6(f), dominates over that in the kinetic energy 
density 6pk = 4>6(f)/a'^. 

This is exactly what happens in the uncoupled models here considered during the 
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Figure 3. Evolution of the DE density contrast, 6^, at the center of a spherical 
matter overdensity for different values of the coupling parameter f3 in RP models. The 
comoving scale of the perturbation is cr = 20h^^ Mpc and the mean matter density 
contrast within the radius i? = cr is (5 = 1 at the present time. 
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Figure 4. Same as Figure [3] but using the SUGRA potential. 
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Figure 5. Top panel: evolution of the mean density contrast, Sdm.b of DM (solid line) 
and baryons (dotted line) in the uncoupled SUGRA case. Bottom panel: same as top 
panel but in the coupled case. Notice how the coupling introduces a bias between DM 
and baryons causing fluctuations to grow at different rates. Plots are given for the 
comoving scale a = 20h~^ Mpc. The mean matter density contrast within the radius 
i? = cr is (5 = 1 at the present time. 



tracking regime in matter era. In fact, in order to have \Spv\ > \Spk\ it must be: 
Spy a pv _ a l - w ^ ^ 

where pk = (j)^ /2a? and py = V are the kinetic and potential background energy densities 
of DE and we made use of the tracker solution oc 5</) oc r^/°'+'^ _ Noticing that in this 
regime we also have w = — 2/(a + 2), relation (|T7|) is always satisfied requiring a > 0. 

The above arguments, however, show that in general a local excess of matter is 
necessary but not sufficient to assure the formation of a corresponding DE underdensity, 
the mechanism being related to the behavior of the background energy density of DE, 
e.g. in our specific case, through the relation: 
2 _ 

Pv > -pk (18) 

a 

Therefore, some differences in the evolution of 5^ arise when DM-DE coupling is 
considered, mainly due to the presence of the 0MDE. As long as the 0MDE holds, the 
kinetic energy of DE dominates over its potential energy invalidating (ITS]) and yielding 
a positive 5^. This is shown in Figures [3] and H] which plot the time evolution of 5^ for 
different values of the coupling parameter j3 in RP and SUGRA cases. 
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Figure 6. Density profiles of DM, baryons (top) and DE (bottom) at the present 
time in tlie uncoupled (left) and coupled (right) RP model. Plots are given for the 
comoving scale a — 20h~^ Mpc. The mean matter density contrast within the radius 
i? = cr is (5 = 1 at the present time. 




Figure 7. Same as Figure [6] but using the SUGRA potential. 



Formation of DE overdensities during (pMDE is understood as follows. Let Sw be 
the first order correction to the state parameter w due to the perturbation in the scalar 
field: 

5w = — (SP^ - w5p^) (19) 
P<t> 

so that w + 6w is the perturbed state parameter inside the fluctuation. Noticing that if 
Pv is negligible with respect to pk, yielding w ~ 1, and 6pk > 0, 6pv < as previously 
observed, we will also have that pk + Spk » pv + Spv- Therefore, to + 5w ~ w ~ 1 (or 
6w ~ 0) and from (fT9|) it follows \Spk\ >> \Spv\- 
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Figure 8. Evolution of the DE density contrast, S^, at the center of a spherical matter 
overdensity in coupled RP models when different scales a are considered. Perturbations 
on shorter scales are suppressed. 
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Figure 9. Same as Figure [5] but using SUGRA potential. 
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After 0MDE, 6^ falls off to negative values resuming the same qualitative behavior 
as in the uncoupled case. Notice that fluctuation amplitudes decrease when increasing 
the coupling strength. In fact, it shortens the period between the end of 0MDE and DE 
dominance so that less time will be left for the perturbation growth. 

As soon as DE becomes the dominant component, starts to decrease to zero as 
the Universe approaches the de Sitter attractor. 

Figure [5] shows the mean density contrast of DM and baryons, 6dm,b, within the 
radius r = cr, as a function of the scale factor. Results are given for the SUGRA model, 
only slightly differences occur when considering different cases. As expected, growth 
of matter perturbations is suppressed once DE dominates the cosmic expansion. Also, 
notice how coupling introduces a bias between baryons and DM causing fluctuations to 
grow at different rates. 

Mean density contrasts, 6dm,b, are easily obtained from ( |T6i) as the assumed initial 
gaussian shape of matter perturbations is kept during their evolution. This is confirmed 
by numerical results and shown in the top panels of Figures (EI) and which display 
the density profiles of DM and baryons for different models at the present time. Density 
profiles of DE are then shown in the bottom panels of the same figures. 

Figures ([8]) and (Q compare the growth of DE perturbations, when different scales 
are considered showing how perturbations on shorter scales are suppressed. 

So far, we have considered the behavior of DE perturbations in the presence 
of matter overdensities. We now look at what happens in the presence of matter 
underdensities or voids. It is not surprising to expect that DE perturbations behave in 
an opposite fashion given that local voids increase the local value of H and consequentely 
6(j), 6(f) < 0. This is confirmed by the Figures (fTOl) . (fTTI) and ( fT2l) which display the time 
evolution of 6^ for matter density contrasts and scales typical of voids and supervoids (see 
next section). Although different models are considered, our results are quite general. 
We find the largest DE inhomogeneities {\6^\ ~ C(10~^ — 10"^)) to be associated with 
objects on very large scale (~ 300/i~^ Mpc) which existence was recently postulated in 
[25] and [26]. On typical supervoid scales we obtain l^^l ~ O{10~^ — 10~^) while we find 
even smaller 15^1 ~ O{10~^ — 10~^) in the presence of voids and supercluster. 

4. Summary and conclusions 

Unlike cosmological constant models (ACDM), models of DE which have a dynamical 
nature yield a varying state parameter w{a). Although current limits on w are consistent 
with a cosmological constant [w = —1 ± 0.1, [2U\), detecting either w 7^ — 1 or its time 
variation {dw/da 7^ 0), would provide a crucial support for dynamical DE. Nevertheless, 
many models of dynamical DE predict no substantial deviation from w = —1 in the late 
time evolution providing a background cosmology very closed to that of ACDM. 

Anyway, dynamical DE stops to mimic a cosmological constant when one deals with 
their clustering properties as dynamical DE is expected not to be perfectly homogeneous. 
Clustering properties of dynamical DE in the vicinity of matter inhomogeneities was 
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Figure 10. Evolution of tlie DE density contrast, (5^, at the center of a spherical 
matter underdensity in RP (top) and SUGRA (bottom) models. Plots are given for 
the comoving scale <j — 2Qh~^ Mpc. The mean matter density contrast within the 
radius R = a is S = —1 at the present time. Such features are typical of matter voids. 
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Figure 11. Same as Figure [TOl but for a — 62h~^ and 5 — —0.75. Such values has 
been estimated for the Bootes supervoid (see text). 
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Figure 12. Same as Figure [TOl but for a ~ 3Q0h^^ and S = —0.3. The existence of 
matter voids with these properties has been recently postulated to explain the observed 
large angle anomalies of the CMB (see text). 

recently studied by Dutta & Maor [9] and Mota et al [10] using two different approaches. 
While the formers face the problem by numerical methods evolving the linearized 
perturbations equations, the latters use an analitycal approach which permitted them 
to extend the analysis even to the nonlinear regime. 

One of the most striking results of the above works is that DE tends to form 
voids in correspondence of linear matter overdensities. However, discrepancies between 
quantitative results are found. When 6 m ~ ^(1); Dutta & Maor [9] find significant void 
amplitudes, ~ (9(10^^), which could be relevant to observations, e.g. on supercluster 
scales. On the other hand, Mota et al [lO] report l^^l ~ (9(10~^) . 

In both the works, DE is modeled as a self-interacting scalar field minimally coupled 
to matter and gravity. 

In this paper we have extended the analysis to a wider class of DE scalar fields 
admitting tracker solutions also allowing for a possible coupling between DM and DE. By 
using the same numerical approach as Dutta & Maor |i9j , we have studied the clustering 
properties of DE in presence of linear matter inhomogeneities on supercluster, void and 
supervoid scales. 

Superclusters are the largest known gravitationally bound massive structures with 
typical radii of about 10 — 25h~^ Mpc and mean density contrasts 5 ~ 1 — 15. Typical 
examples are the local supercluster (LSC), of which our galaxy is a part, and the Shapley 
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supercluster (SSC) which is ~ 650 Mlys from us. The LSC has a mean overdensity 
of 5m ~ 2 — 3 over a scale ~ 15h~^ Mpc |2T] while the SSC has been found to 
have an overdensity of 6m ~ 10.3 over a scale of 10.1h~^ Mpc [22]. In addition to 
matter overdensities, Universe contains voids of matter typically having radii similar to 
superclusters or even larger (supervoids). Data from 2dFGRS are consistent with voids 
having average radii of ~ 15/i~^ Mpc and average mean density contrast 6^ ~ —0.93 |23j . 
An example of a supervoid is the Bootes void found to be rougthly spherical with radius 
of ~ 62h~^ Mpc |2l] and a mean density contrast estimeted to be —0.8 < 6m < —0.66. 

Recently, the possibility that extremely large voids might exist with radii of 
~ 100 — 300 Mpc, has been considered. In particular in [25] it is shown that the 
Integrated Sachs Wolfe (ISW) effect due to a void with radius ~ 200 — 300 Mpc and 
6m ^ — 0.3 would be observed as a cold spot in the CMB radiation explaining the 
observed large angle CMB anomalies (see also [26] for a similar conclusion). According 
to our results (and those of Mota et al [10] as well) such extremely large objects would 
correspond to the largest DE overdensities. 

We conclude summarizing our main results concerning the behavior of DE 
perturbations in the vicinity of matter inhomogeneities. In the presence of matter 
overdensities we find that: 

(i) DE tends to form voids if no coupling to DM is present 

(ii) in coupled models, DE overdensities form during 0MDE. After this stage (5</, becomes 
negative resuming the same behavior as in the uncoupled case. 

If matter underdensities are considered, DE perturbations behave in the opposite 
fashion. From our analysis we obtain |5</,| ~ O{10^^ — 10^^) for voids and supercluster, 
\6(j,\ ~ C(10~^ — 10^^) for supervoids and if extremely large voids exist ~ 
(9(10"'^ — 10^^). Our results are consistent with those of Mota et al |10j indicating 
that if DE is described by a scalar field, either uncoupled or coupled with DM, it would 
be almost homogeneous on sub-horizon scales above that of galaxy clusters. 

Further, in general we have: 

(i) within matter inhomogeneities DE can form voids as well as overdensities. However, 
the behavior of 5^ is mainly related to how the background DE energy density 
evolves rather than to whether a matter fluctuation is an overdensity or an 
underdensity 

(ii) DE perturbation growth is sensitive to the scale considered. On shorter scales 
perturbations are suppressed 

(iii) accelerated expansion yields a suppression of DE inhomogeneities, i.e \6^\ decreases 
to zero as the universe approaches the de Sitter attractor. 

We have shown in the previous section that decreases at the increasing of the 
coupling strengh (3 when supercluster and void scales are considered. According to the 
first of the above claims, the fluctuation suppression due to the coupling is then to impute 
to the fact that perturbations grow on horseback of two distinct evolutionary phases of 
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the scalar field, namely 0MDE and usual tracking phase in matter era, changing their 
sign. It would be interesting to ascertain whether such effect also holds when dealing 
with nonlinear matter collapse, e.g. on galaxy cluster scales. If so, very small DE 
perturbations are expected, at most not larger than what found by Mota et al [10] when 
considering nonlinear scales in uncoupled DE models, i.e. 6^ ~ O{10~^). This would be 
a quite unexpected result since it seems natural to believe that infalling matter, when 
coupled to DE, will drag along DE permitting larger DE perturbations. 

Suppression of DE perturbations due to the coupling is however strictly true 
only after 0MDE if larger scales are considered. During 0MDE, DE fluctuations are 
only marginally suppressed on supervoid scales, while they can be as larger as matter 
fluctuations on extremely large scales (~ 200 — 300 Mpc). Figure [T2] indicates 5,^ ~ 6m 
at z 100 and even larger values are expected for higher z. suggesting the possibility 
that DE clustering might be detected through the ISW effect. 

As already observed in [27], an interaction between DM and DE changes both the 
scaling of the DM energy density and the growth rate of matter perturbations affecting 
the time evolution of the metric potentials and, consequentely, the ISW effect. Our 
results point out that, in coupled cosmologies, a further contribution to ISW effect can 
arise during 0MDE from DE perturbations associated with very large voids of matter. 

Further investigations on this last point and the behavior of DE in the presence of 
nonlinear matter inhomogeneities are left to future works. 
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Appendix A. 

In a spacetime described by the metric (j4]): 
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(A.6) 
(A.7) 
(A.8) 
(A.9) 

Here, overdots and primes denote the derivatives with respect to t and the radial 
coordinate r respectively. Only radial motions are considered and terms quadratic in 
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has been neglected as DM and baryons are non-relativistic components. 



According to ([5]) and (Ej) we decompose our variables in an homogeneous part and 
a perturbation. It is then straightforward to obtain the equations for the background: 



37^2 = SttG [pdm + Pfe + P0] 
'4> + m'^ + V' = Cpam 

Pdm + '^T^Pdm = —Cpdm<P 
Pf, + SHpb = 

To linear order, ( ]A.2l) -( lA75l) gives the equations for the metric perturbations: 

X' 2^' 



ip" = 4nG {5 Pdm + 5pb + 



C + 2ff(2C + ^)+^(- 
^ + H (hij + C) + - ^ + - - — - i^''] =^TrG[5pdm + 5pb + l 

2iJ (C + 2^) + ^ f 4 - 4 + - - — - ^"1 = SttG {5pdm + 5pb + 5p^) 



-{C-i,)-2i,' 



(A.ll) 
(A.12) 
(A.13) 
(A. 14) 

(A.15) 

(A.16) 
(A.17) 



PdmVdm + PbVl + [pcf, + P<^) Vl 



while perturbation equations for DE, DM and baryons follow from (lA.6p - (lA.10p : 
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6pb + 3/75p, + (C + 2^) Pb + (^i;r + j Pb = 
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(A.18) 

(A.19) 

(A.20) 

(A.21) 

(A.22) 
(A.23) 
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Combining equations (lA.lSp . (1A.16P and (1A.17P gives: 

(C + 2ij) +2H[c + 2ij) = -47iG {6pdm + Spb + 6p^ + 36 P^) (A. 24) 

It is clear, from the above equations, that x = C + 2"?^ is the only combination which is 
relevant for the evolution equations of DE, DM and baryons perturbations. 

When rewritten in terms of the conformal time r and the variables Xy = 6pi/pi 
and Oi = vl' + 2</r {i = dm, b, (p) equations fjXTT]) - flA:i4l) and ^KM-^KMi gives (El) 
and ([8]) respectively. 
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